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Aims. On the basis of the PPMXL catalogue we perform an all-sky census of the Hyades down to masses of about 0.2 nio in a region 
up to 30 pc from the cluster centre. 

Methods. We use the proper motions from PPMXL in the convergent point method to determine probable kinematic members. From 
2MASS photometry and CMC14 r'-band photometry, we derive empirical colour-absolute magnitude diagrams and, finally, determine 
photometric membership for all kinematic candidates. 

Results. This is the first deep (r' < 17) all-sky survey of the Hyades allowing a full three-dimensional analysis of the cluster. The 
survey is complete down to at least Mk, = 7.3 or 0.25 hiq. We find 724 stellar systems co-moving with the bulk Hyades space velocity, 
which represent a total mass of 435 mg. The tidal radius is about 9 pc, and 275 mg (364 systems) are gravitationally bound. This is the 
cluster proper. Its mass density profile is perfectly fitted by a Plummer model with a central density of 2.21 mgpc"^^ and a core radius 
of r^o = 3.10 pc, while the half-mass radius is rn = 4.1 pc. There are another 100 mg in a volume between one and two tidal radii 
(halo), and another 60 mg up to a distance of 30 pc from the centre. Strong mass segregation is inherent in the cluster. The present-day 
luminosity and mass functions are noticeably different in various parts of the cluster (core, corona, halo, and co-movers). They are 
strongly evolved compared to presently favoured initial mass functions. The analysis of the velocity dispersion of the cluster shows 
that about 20% of its members must be binaries. As a by-product, we find that presently available theoretical isochrones are not able 
to adequately describe the near-infrared colour-absolute magnitude relation for those cluster stars that are less massive than about 0.6 
mg. 

Key words, open clusters and associations: individual: Hyades; Stars: luminosity function, mass function; Hertzsprung-Russell and 
C-M diagrams 



1. Introduction 

Without doubt, the Hyades are one of the best-studied open 
clusters in our Galaxy. Its proximity to the Sun has made it 
an interesting target for centuries, so it is impossible here to 
give appropriate credit to all the astronomers who worked on 
the Hyades. With no claim to c o mpleteness, we m e ntion here 
the stu d ies bvlvaii BuerenI (Il952l).lvan Altenal (ll969h.lPels etalj 
(Il975h . iHansonl (Il975h lGunnetal.1 (Il988l) . and iReidI (Il992h . 



Hogan et a 0112008!) have surveyed some 275 deg2 of 



Recently, 

the Hyades based on a combination of 2MASS and UKIDSS 
observations. 

Loosely speaking, we distinguish between pencil beam stud- 
ies and those relying on an all-sky survey. All the above- 
mentioned studies are pencil beam studies since they only cover 
a limit ed field-of-view, ev en if it is as large as 275 deg^ in the 
case of lHogan et alj ( l2008h . On the sphere, this only corresponds 
to a distance of 7.3 pc from the Hyades centre. Pencil beam stud- 
ies are based on dedicated observations and may comprise a vari- 
ety of astrophysical parameters. They can only indirectly reveal 
the full three-dimensional structure of the cluster. 



Send offprint requests to: S. Roser 



An all-sky survey allows studying the cluster in its full three- 
dimensional extent, but one usually relies on the few avail- 
able data entries in the survey, e.g. proper motions or photom- 
etry in a few bands. Of these studies we mention the work by 
Schwan (1991) using the old PPM and FK5 catalogues and the 
papers bas ed on the Hipparcos obs ervations by Perrvm an et al.l 
(Il998h and lde Bruiine et al.l (l200lh . These all-sky studies were 
restricted to stars brighter than about V » 11 to 12 (Hipparcos) 
without being complete to this magnitude. 

As far as membership of individual stars in the Hyades 
cluster is concerned, the results of the studies above were col- 
lected in the Prosser & StaufFer data base (presently available 
from J. StaufFer, priv. comm.) some fifteen years ago. Prosser & 
Stauffer's data base was assumed to be complete down to 0.1 
m© (Bouvier et al., 2008). Much effort has since been dedicated 
to the search for very low-mass (below 0. 1 m©) H yades candi- 
dates. These attempts have been partly unsuccessful (iGizis et al.L 
1999), or resul ted in the confirmation of membersh i p of on ly one 
M8.5 dwarf (Reid & Hawlevl lT999l) . lHogan et"an ( 120081) report 
the detection of 12 L-dwarfs in a 275 deg2 field of view, and 
iBouvier et al.l (l2008h describe the detection of two brown dwarfs 
in an area of 16 deg^ around the Hyades centre. 
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Using PPMXL ( iRoeser. Demleitner & SchilbachLl2()Toh we 
can, for the first time, extend previous studies by going as deep 
in magnitude as the pencil beam studies and still performing 
an all-sky survey. We supply a complete sample of candidates 
within 30 pc from the centre, determine individual distances to 
the stars, and so resolve the full three-dimensional spatial struc- 
ture of the Hyades, and its immediate neighbourhood, from the 
most massive stars down to almost 0. 1 m© by an application of 
the convergent method to the proper motions in PPMXL. 

The paper has the following structure: in the next section 
we briefly introduce the observational dataset from PPMXL we 
are using. In section|3]we describe the convergent point method. 
Then, in section |4] we define our criteria for membership deter- 
mination and discuss the amount of possible contamination of 
the sample. We compare the empirically derived colour-absolute 
magnitude diagram of the Hyades with theoretical isochrones in 
section 14.41 In section |5] we reveal the spatial structure of the 
Hyades and derive the shape parameters of the cluster Section 
|6]is devoted to the present-day luminosity and mass functions. 
Then we discuss the internal velocity dispersion in section |7] 
and, finally we summarise our results in section |8] 



2. Observations 

As primary data of observations w e use the PPMXL catalogue 
(iRoeser. Demleitner & Schilbachi l2010h . PPMXL contains 
positions and proper motions on the ICRS for some 900 million 
stars (from the brightest stars down to about V 20) and covers 
the complete sky. PPMXL also gives low-accuracy photometry 
from USNO-BLO (Monet et al., 2003), a subset of 410 million 
stars contains 2MASS (iSkrutskie et al.l i2006l) photometry. The 
typical individual mean errors of the proper motions range 
from better than 2 mas/y for the brightest stars with Tycho-2 
(lH0g et al.L l2000h observations to more than 10 mas/y for the 
faintest stars in the region south of -30° declination. For the 2.5 
million brightest stars Tycho-2 B, V photometry is available. 

For this work we cr oss-ma tched PPMXL 
with UCAC3 (Zachari asetalL 1201 Oh and CMC14 
( [Copenhagen Univ. Obs. et al. L I2006h . We combined the 
positions and proper motions of PPMXL with the positions 
from CMC14 and from UCAC3. The latter had to be recon- 
structed because the original positions of UCAC3 are not 
published. With these data we performed a weighted least- 
squares adjustment to derive new, improved, mean positions 
and proper motions. In the following we refer to this subset of 
PPMXL as the Carlsberg-UCAC (CU) subset. It turned out that 
the most important feature in the CU is the accurate photometry 
in at least one optical band from CMC 14, the r'-magnitudes in 
the SDSS system. CMC14 observed the sky in the declination 
range from -30 deg to about +52 deg, and it is 95% complete 
down to r' - 16.8 and 80% to 17.0. Then completeness drops 
rapidly (when compared to SDSS). Its limiting magnitude 
is r' = 17.8 (ICopenhagen Univ. Obs. et al. L l2006l) . The CU 
contains some 140 million stars, 90 million of which have 
r'-magnitudes from CMC14. 

The astrometric information for the 120 000 Hipparcos 
stars is ta ken from the new reduction of the Hipparcos data by 
Ivan Leeu wen (2007) instead of the data contained in PPMXL, 
the cross-matches with 2MASS are kept as are the B, V magni- 
tudes from Tycho-2 of the Hipparcos stars. 



3. Convergent point metliod 

For a nearby open cluster like the Hyades the convergent point 
method is a suitable tool to determine membership when only 
proper motions and (for a representative subset) radial veloci- 
ties ar e avail able. The method is textbook knowledge (see, e.g. 
ISmartlll968h, an d an e xcellent recent description can be found 
m van Leeuwenl(l2009h . Hence, we will not repeat a description 
here. 

In this paper, we did not apply the convergent point method 
from scratch, but used the six-dimensional phase space parame- 
ters fo r the Hyades cluster centre from Table 7 in Ivan LeeuwenI 
(I2009h . We use the galactic rectangular coordinate system X, Y, Z 
with origin in the Sun, and axes pointing to the Galactic Centre 
(X), to the direction of galactic rotation (F), and the North 
Galactic Pole (Z). In galactic coordinates, we adopt [Xc, yc, zd 
- [-43.1, 0.7, -17.3] pc for the position and [u^, Vc, wd = [-41.1, 
-19.2, -1.4] kms ' for the mean motion of the cluster centre. In 
other words, the stars most accurately measured, Hipparcos stars 
from the new reduction by van Leeuwen (2007,) , determine the 
position and the motion of the Hyades cluster. The n the direction 
of the convergent point is also given, see Eq. 7 in Ivan LeeuwenI 
(I2OO9I) . 

Given the phase space parameters from van Leeuwen '(' 2009I) . 
we search for all stars in the sky (in this case the PPMXL cata- 
logue, resp. the CU subset) which have proper motions consis- 
tent with the given space velocity Vf - [Uc, Vc, Wc] of the cluster 
centre. For each line-of-sight (a, 6) we can then express the vec- 
tor of the space velocity by one component Vi parallel to the 
line-of-sight, the radial velocity, and one perpendicular to it Vt, 
the tangential velocity (in the tangential plane). It is easy to see 
that, for a Hyades member, the component of Vc parallel to the 
line-of-sight only depends on (a, 6), and not on the distance of 
a star from the Sun. Hence, the same is true for the component 
of Vf in the tangential plane. The tangential velocity can be split 
into one component vy in the direction to the convergent point, 
and a component Vj^ perpendicular to it. The expectation value 
of the latter is zero. 

The PPMXL catalogue gives for each star the positions (a, 6) 
and the proper motions (fj.a,l-is), i e- four of the six components 
of the phase space coordinates, the radial velocity Vr and the dis- 
tance d from the Sun being the remaining two. A necessary con- 
dition for a star to be a kinematic member of the Hyades cluster 
is: 



|Vc - v| < e. 



(1) 



where v = [u,v,w] is the space velocity of a candidate and e is a 
suitable bound. As long as the radial velocity Vr is not measured, 
we adopt as Vi the projection of the cluster motion Vc to the line- 
of-sight. 

At each position {a, 6) the proper motions {iJa,f^s) can be 
rotated into the direction of vy and Vj^. We call these motion 
components /iy and ju±. The remaining parameter, distance d or 
parallax w, is given by 



fx 



V|| 



(2) 



where / - 4.74kms"'kpc"' and m is the so-called secular par- 
allax. As the method fixes v, as explained above and vj from 
Eq. 2, the necessary condition (Eq. 1) for kinematic membership 
shrinks to |vx| < e. Since the expectation value of v^ is zero, 
< (v^) >2 is a measure for the one-dimensional velocity disper- 
sion of the cluster. 



2 



S. Roser et al.: A deep all-sky census of the Hyades 



The convergent point method predicts the radial velocity v, 
and the distance d from the Sun for a cluster candidate. A fi- 
nal confirmation as a member can be obtained by directly mea- 
suring Vi and m. Other than from trigonometric parallaxes, the 
distance d from the Sun can be confirmed also by determining 
photometric parallaxes. The CU subset of PPMXL catalogue 
gives 2MASS J, H, /Tj-photometry plus r'-magnitudes. For the 
2.5 million brightest stars in PPMXL B and V magnitudes from 
Tycho-2 are available. These photometric data are used to con- 
firm or reject photometrically the parallaxes from the convergent 
point method. The radial velocities v, can only be confirmed by 
direct measurements. 



4. Membership determination 



4.1. Kinematic membership 

The Hyades have a tidal radius of about 10 pc (see, e.g. 
|Perrvmanetanil998h and an average velocity dispersion in one 
dimension of 0.23 kms"' ([Ounn et al., 1988). In this paper we do 
not concentrate only on the gravitationally bound cluster itself, 
it is also interesting to investigate the situation in the immediate 
surroundings. There we expect to reveal former Hyades mem- 
bers which already left the cluster because of dynamical evolu- 
tion. 

Therefore, we allow a bound e - 4kms ' for the velocity 
component Vj.. Additionally, we have to take into account that 
V|| goes to zero when (a, 5) approaches the convergent point. 
This led us to set an upper bound 77 for the tangent of the an- 
gle between the proper motion vector and the direction to the 
convergent point 

^ = ^ < 77. 

At the cluster centre rj - 116 corresponds to |vj^| - 4kms"', 
so we set this condition to hold everywhere in the cluster 

We also specified a third bound by requiring that the distance 
rc of a candidate from the cluster centre should be less than 30 
pc which is sufficiently large compared to the tidal radius. 

These kinematic selection criteria are fulfilled by 15757 stars 
out of the 140 million contained in the CU subset. They are sub- 
ject to the photometric selection described in the next section. At 
this stage we excluded all candidate white dwarfs from further 
consideration. White dwarfs in the Hyades are at the limit of be- 
ing observed in 2MASS, and if so, the photometric accuracy is 
low. That prevents us from using photometric distances of white 
dwarfs from the CU to check their predicted secular parallaxes. 
We will investigate the white dwarfs associated with the Hyades 
in a separate study. 



hence reddenning has been neglected in all future considera- 
tions. 



4.2. Photometric membership 

In the following we will check if the kinematically selected 
stars populate allowed loci in the colour-absolute-magnitude di- 
agrams. As the candidates occupy a large magnitude range, there 
is no unique CMD where all the stars have sufficiently accurate 
photometry. Stars brighter tha n V - 10 have good quality in 
B and V mainly from Tycho-2 (lH0g et al.L 12000 ). or ASCC-2.5 
(lKharchenkoLl200lh . On the other hand, stars fainter than about 
Ks = 5 have good quality photometry from 2MASS and CMC14. 

Fortunately, the interstellar reddenning towards the Hyades 
is very low . In a critical r eview on reddening determinations of 
the Hvades lT^vlo^ (l2006l) concludes that E{B - V) < 1 .0 mmag. 




Fig. 1. The Mk, vs. r' - Ks diagram of the kinematic candidates 
(small dots). The vast majority of the stars have kinematically 
predicted parallaxes inconsistent with their photometric paral- 
laxes. The candidates shown as thick dots are retained for fur- 
ther analysis. Stars brighter than Mk, ~ 3.5 are not contained in 
CMC 14. 



For the absolutely brightest stars we rely on the classical 
My VS. B - y diagram for cross-checking the kinematically pre- 
dicted parallax m with the photometric parallax. In fact, almost 
all the bright Hyades members have been obser ved by Hipparcos 
and their membership is settled (see, e.g. P errvman et al.l IT998I: 
[de Bruijne et al., 2001; van Leeuwen, 2009'). 

For stars with masses less than about 0.6 m© and the age and 
metallicity of the Hyades there are discrepancies between dif- 
ferent theoretical isochrones, especially in the NIR. We discuss 
this in more detail in Sect. 14.41 As a consequence of this, we 
have no choice but to select the Hyades candidates from an em- 
pirically derived colour-magnitude diagram obtained from the 
distribution of the kinematically selected candidates. In the Mk, 
vs. J - Ks diagram the selection of low-mass stars (< 0.6 m©) is 
hampered by the fact that the colour J - K, is nearly constant. 
In a first step, we therefore used the M^^ vs. r' - (FiglT]l and 
the Mk^ vs. J -Ks (Fig.|2]i diagrams to discard the large majority 
of the 15757 kinematic candidates whose photometric parallaxes 
were inconsistent with the predicted secular ones. 

In order to work out more clearly the empirical Hyades 
main sequence, we selected, in FigslU and |2] only stars with 
Iv^l < Ikms"', i.e. the k,9QQQ most probable kinematic can- 
didates. In FiglT] we show the Mk, vs. r' - Ks diagram of all 
those 9000 stars (small dots) that have r' magnitudes, i.e. that 
are contained in CMC14. The stars brighter than Mk, ~ 3.5 
are missing. They are not contained in CMC14, but their pho- 
tometric membership is determined from Tycho-2 B, V photom- 
etry. The Hyades main sequence stands out prominently, and the 
candidates are provisionally marked as thicker dots. Stars below 
the Hyades sequence appear as sub-luminous stars and must be 
much farther away than Hyades members. In other words, their 
parallax m predicted from the kinematic selection is inconsistent 
with the photometric parallax, and hence these stars cannot be- 
long to the Hyades. In the Mk, vs. J - Ks diagram for the same 
stars (Fig|2]i we find a considerable density of stars at colours 
J - Ks between 0.8 and 1. We marked the candidate Hyades 
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from Fig. [Has thick dots also in Fig.|2] It is clearly visible that 
NIR photometry alone is not sufficient to properly select faint 
stars with Mk^ > 5: from Fig. [T] we find no Hyades candidates 
with Mk, > 9 although Fig|2] suggests plenty of them. Instead, 
all the candidates with Mk^ > 9 have r' - colours which are 
to blue, and this means that their predicted secular parallaxes are 
too large. We conclude that stars with Mfc^ > 9 in the Hyades are 
below the limiting magnitude of CMC 14. 

Even if we consider the finite width of the Hyades sequence 
caused by the presence of binaries, we still find a number of stars 
above the sequence in both CMDs (FiglUand Fig.|2]l. In their ma- 
jority they are field giants or very seldomly foreground dwarfs. 
The former are discarded via the two-colour, J - H \s. H - Kg, 
diagram. In the Mk^. vs. J - Kg diagram we note a few white 
dwarfs near J - ~ and 11 < Mk^ < 14, but investigating 
them needs a different approach. 




-05 00 05 10 15 

J-K s 

Fig. 2. The Mk, vs. J-Ks diagram of the candidates. Stars which 
were marked as thick dots in Fig. [T]are also marked as thick dots 
here. 

For the final selection we determined the lower bound of an 
empirical cluster sequence stepwise over r' - and J - Ks, and 
selected all stars between this lower bound and the upper limit 
given by binaries. Out of the 15757 kinematically selected candi- 
dates, we could confirm the secular parallax via the photometric 
parallax for 775 stars. 

4.3. Final membership selection 

The sample of 775 kinematic and photometric candidates was 
cross-matched with the probable Hyades members from the 
Prosser & Staufifer database. Out of their 536 candidates, 407 
were identified in our sample. Further, we found 6 stars which 
belong to the Hyades with certainty but have been missing in 
our list due to erroneous proper motions in PPMXL (checked 
via Vizier and Aladin). These stars were added to our sample to 
give a total of 781. The membership of the majority of the re- 
maining stars from the Prosser & StaufFer database has not been 
confirmed since they did not pass either the kinematic or photo- 
metric criteria. A few stars from Prosser & StaufFer are too faint 
to be observed in CMC14, so we could not check their member- 
ship in a proper way. 

4.3.1. Contamination estimate 

Although the kinematic and photometric criteria provide a pow- 
erful tool for selecting cluster candidates and enable us to isolate 
781 candidates out of 140 million on the whole sky, we cannot 



exclude a certain contamination of our sample by field stars. The 
contamination can be estimated empirically by directly compar- 
ing the predicted parallaxes and radial velocities with the cor- 
responding measurements as far as these are available. Another 
approach is to compute the probability of contamination from a 
kinematic model of the Galaxy, or, alternatively, by considering 
the observed velocity dispersion of field stars in the solar vicin- 
ity. In the following we apply these methods to estimate probable 
contamination of our Hyades sample. 

For about 300 stars from the candidate list we could find 
Hipparcos parallaxes and/or radial velocities in Vizier/Simbad. 
We call them "control stars" in the following. Whereas for the 
majority of the control stars the predicted and measured param- 
eters coincide well within the standard errors, 26 of them have 
observed parallaxes and/or radial velocities that differ signifi- 
cantly from the computed ones. We analysed their space co- 
ordinates predicted by the convergent point method and found 
that all of them are relatively far away from the cluster centre. 
Indeed, the membership was not confirmed for all control stars 
with Iz-Zcl > 20 pc (1 1 stars). The remaining 15 "wrong" control 
stars are located at distances r^ from the cluster centre larger than 
9 pc. From their distribution we estimate a 7.5% contamination 
at 9 pc < rc < 18 pc, 30% at 18 pc < r^ < 30 pc and 100% when- 
ever |z - Zcl > 20 pc. Following this finding, we exclude all stars 
from the Hyades candidate sample which are located more than 
20 pc away from the cluster centre in z-direction. This leaves 
us with a sample of 724 candidates, and the formal application 
of the empirical rule found above predicts some 65 field stars 
among these 724 Hyades candidates. 
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Fig. 3. The luminosity functions (LF) of the Hyades candidates 
and of field stars in the solar vicinity, a: Solid lines and symbols 
are the LFs of the Hyades at different distances from the cluster 
centre rc, black filled circles for rc < 9 pc, red open circles for 
9 pc < rc < 18 pc, blue crosses for 18 pc < rc < 30 pc. The 
LF of nearby stars from the CNS4 (JahreiB, H. , priv. comm.) is 
shown by the black dashed line whereas the black dotted line is 
the fraction of field stars with space velocities comparable to the 
velocity of the Hyades. b: The red open circles and blue crosses 
show the adopted contamination of the Hyades sample of 7.5% 
at 9 pc < rc < 18 pc and 30% at 18 pc < rc < 30 pc, respectively. 
The black dotted line is the same as in a. 

Since the control stars are relatively bright stars (see Fig.|4| 
that usually have astrometric and photometric data of better ac- 
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curacy, we must consider the possibility that the contamina- 
tion rate may be somewhat higher for the faintest stars of our 
sample. Therefore, we checked whether our empirical estimate 
from above coincides with the statistical predictions of galac- 
tic models within the complete magnitude range covered by 
our sample. We used the kinem atic model of the Galaxy by 
iKharchenko & Schilbachl(ll996h to estimate the number of field 
stars within a circle of a radius 5 deg from the cluster centre and 
with proper motions compatible with Hyades membership in this 
region of the sky. The model predicts 5.06 such field stars in an 
area of 78.54 deg^ in the direction of the Hyades. That suits per- 
fectly the number of 4.95 stars we obtain if we use the empirical 
rule of the contamination estimates from above. 

Further, we computed the probability for field stars in the 
local solar neighbourhood to have, by chance, the same space 
velocities as the Hyades members. JahreiB (priv. comm.) has re- 
cently determined the K-band luminosity function of the stars in 
the solar neighbourhood from the CNS4, the new 4th Catalog 
of Nearby Stars (in preparation) from which we infer the vol- 
ume density of the local sample. The CNS4 luminosity function 
is shown in Fig. [3^ together with the luminosity functions of 
the 724 Hyades candidates at different distances from the clus- 
ter centre. The random contamination of our Hyades sample 
is now determined by the percentage of stars having the same 
space velocity [u, v, w] as the Hyades members with respect to 
the l ocal standard of rest (L SR). Corrected for the solar motion 
from iPiskunov et aP (|2P06), the LSR velocity of the Hyades is 
[-31.8, -7.1, -1-6.1]. Assuming that the velocity distribution of 
field stars follows a three-dimensional normal distribution with 
standard deviations given by the velocity dispersions, we com- 
pute the number of field stars with H yades motion in the range 
< Mk^ < 8. We used the data from ioehnen & BinnevI d 19981) 
who determined the velocity dispersions of field stars of differ- 
ent colours in the solar neighbourhood from Hipparcos measure- 
ments. Taking into account that only two velocity components 
are observed, we find that the percentage of random field stars 
with Hyades motion is less than 0.6 x 10"^ of the local density, 
and this is plotted as the dotted line at the bottom of Fig. [3^. 

Since the spatial density of the Hyades candidates with r^ < 
9 pc is, on average, by two orders of magnitude higher than that 
expected from the local density of field stars, we conclude that 
contamination is negligible for this volume. With increasing dis- 
tance from the cluster centre, the density of the Hyades stars 
is decreasing, and the probability of contamination becomes 
higher For comparison we show in Figj3j) the contribution of 
field stars (the same as in Fig. |3^) together with the contamina- 
tion expected in outer regions, like what we obtained empirically 
a few paragraphs above. For 9 pc < r^ < 18 pc, the contamina- 
tion by field stars is always less than the estimated 7.5%, except 
for the faintest stars with Mk, > 8 where our catalogue becomes 
incomplete. Also in the region 18 pc < rc < 30 pc, the contam- 
ination of 30% from our empirical finding coincides well with 
the statistical estimates in this magnitude range. 

To summarise: in the following, our study is based on a sam- 
ple of 724 Hyades candidates. The contamination by field stars 
is assumed to be dependent on the distance rc from the cluster 
centre: a negligible contamination at rc < 9 pc, 7.5% at 9 pc 
< rc < 18 pc, and 30% at 18 pc< rc < 30 pc. This is taken into 
account when we discuss the density distribution in the Hyades 
cluster, its luminosity and mass functions. 



4.3.2. The final sample 

The 724 Hyades candidates are listed in Table [T] organised as 
follows: in the first column we give a running number The next 
six columns are the right ascension and declination for equinox 
and epoch J2000.0, the corresponding proper motions, and the 
mean errors of the proper motions. The following 12 columns (8 
to 19) contain photometric data, i.e. magnitudes and their cor- 
responding mean errors: Johnson B and V, mainly from ASCC- 
2.5, /, H and K, from 2MASS, r' from CMC14. All the follow- 
ing columns are derived quantities. Columns 20 and 21 are the 
secular parallax m and its mean error, column 22 gives the pre- 
dicted distance from the cluster centre, column 23 is the pre- 
dicted radial velocity, columns 24 to 27 are vy and v^ with cor- 
responding mean errors. In the last column 28 we give the mass 
of the stars as derived in section |5]below. This table is available 
from the CDS, Strasbourg, France. 



Table 1. Description of the table containing the data of 724 
Hyades candidates. This table is available from the CDS, 
Strasbourg, France. 





Label 


Units 


Explanations 


1 


name 




Running number 


2 


RAdeg 


deg 


Right Ascension J2000.0, epoch 2000.0 


3 


DEdeg 


deg 


Declination J2000.0, epoch 2000.0 


4 


pmRA 


mas/yr Proper motion in RA-cos(DE) 


5 


pmDE 


mas/yr Proper motion in DE 


6 


SpmRA 


mas/yr Mean error of pmRA-cos(DE) 


7 




mas/yr Mean error of pmDE 


8 


B 


mag 


Johnson B magnitude from ASCC-2.5 


9 


eg 


mag 


Standard error of B magnitude 


10 


V 


mag 


Johnson V magnitude from ASCC-2.5 


11 


ey 


mag 


Standard error of V magnitude 


12 


J 


mag 


J magnitude from 2MASS 


13 


ej 


mag 


Standard error of ] magnitude 


14 


H 


mag 


H magnitude from 2MASS 


15 


en 


mag 


Standard error of H magnitude 


16 


Ks 


mag 


K, magnitude from 2MASS 


17 




mag 


Standard error of magnitude 


18 


r' 


mag 


SDSS r' magnitude from CAMC14 


19 




mag 


Standard error of r' magnitude 


20 


m 


mas 


Secular parallax 


21 




mas 


Mean error of Secular parallax 


22 


Tc 


pc 


Distance from the cluster centre 


23 


Vr 


kms"' 


Predicted radial velocity 


24 


V|| 


kms"' 


Tangential velocity in direction 








to the convergent point 


25 


ev|| 


kms"' 


Mean error of vy 


26 


v± 


kms"' 


Tangential velocity perpendicular to the direction 








to the convergent point 


27 




kms"' 


Mean error of v^ 


28 


m 


mo 


mass of the star (from M/L relation, 



see Sect.O 



The vs. J - Ks colour-magnitude diagram of the final 
sample is shown in Fig. |4] Since the CU subset of PPMXL is 
95% complete down to r' = 16.8, we find that our sample starts 
getting incomplete between Mk, - 7 and Mk, = 8. A subset of 
stars in Fig. |4]is shown as green dots. For these stars we were 
able to confirm their predicted secular parallaxes and/or radial 
velocities via measured trigonometric parallaxes and/or radial 
velocities. 
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Fig. 4. The Mk, vs. J-Ks diagram of the 724 Hyades candidates. 
For Mk^ < 2, the J - Kg colours from 2MASS are of low preci- 
sion, because the stars are too bright in the J-band. Stars shown 
as green dots are those for which we were able to confirm their 
membership by comparing with measured trigonometric paral- 
laxes and/or radial velocities. 



4.4. Comparison with tiieoretical isociirones 

The empirical selection of kinematic members was chosen 
as present-day theoretical isochrones fail to describe the main 
sequence of the Hyades. We compared the empirical colour- 
magnitude diagram M^^ vs. / - K, with 3 different theoretical 
isochrones. In the literature, we found estimate s of the metallic- 
ity of the Hyades, e.g. as [Fe/H ] - 0.14 ± 0.05 dPerrvman et al.L 
[T998h . [ Fe/HI ^0.144 + 0.013 (lGrenonil200l and [Fe/H] =0.13 
± 0.01 (iPaulson e t al.'. '2003'). The age of the Hyades is given 
as 625 ± 50 Myr (Ferryman et al., 1998) or 648 ± 45, the latter 
based on cooling white dwarfs dPe Gennaro et al.L 12009). From 
models taking i nto ac count stellar rotation of the turn-off stars 
iLebreton etall (l200ll) estimate a range of 500 to 650 Myr for 
the Hyades age. 

In Fig. |5] we compare the empirical colour-magnitude 
diagram with available theoretical isochrones for ages and 
metallicities as close to the Hyades valu e s as p ossible. The 
isochrones ai-e taken from iMarigo et"aD (120081) (PADOVA, 
650 Myr, [Fe/H] = 0.14, filled red triangles), 'Pott er et all 
(2008) (DAR TMOUTH, 600 Myr , [Fe/H] = 0.21, open blue 
circles ) and iBaraffe et all (Il998h (BCAH, 650 Myr, [Fe/H] 
= 0.0, open magenta triangles). In the latter case we used 
BCAH98_models.l. The BCAH models are insensitiv e to the 
L_mix parameter below 0.6 m© (iBaraffe et al.L Il998h . When 
isochrones were not given in 2MASS magnitude bands, we con- 
verted them using the formulae given in the 2MASS webpage 
www . ipac . caltech.edu/2mass/. 

We find that the DARTMOUTH isochrones describe the Mk^ 
vs. J - Ks of the Hyades perfectly in the regime brighter than 
Mk^ - 5.0, but fail to describe the situation for fainter stars. This 
is, Hyades stars with masses less than 0.65 m© have colours red- 
der in y - K, than predicted by DARTMOUTH for stai's with 
corresponding age and metallicity. The less massive the stars 
are, the greater the discrepancy. The PADOVA isochrones also 
coincide well with observations for J - Ks < 0.7 (0.7 m©). For 



Fig. 5. Comparison of the empirical Mk, vs. J-Ks diagram of the 
Hyades candidates (masses less than about 1 m©) with different 
theoretical isochrones. The PADOVA isochrones are shown by 
filled red triangles (curve most to the left), the DARTMOUTH 
isochrones blue open circles, and the BCAH isochrones by the 
magenta open triangles (curve most to the right). 



stars with masses less than 0.7 m© the theoretical J - is about 
0. 1 mag too small. 

Finally, a comparison has been made with the BCAH 
isoschrones. They cover a mass range from 0.06 to 1.1 m©. For 
masses larger than 0.8 m© (J - K, = 0.6), BCAH gives Mk,- 
magnitudes which are slightly too bright compared to PADOVA 
and DARTMOUTH. For 4.5 < Mk, < 6.5 the theoretical J-K^- 
colours are slightly too blue compared to observations. But for 
fainter stars the isochrones satisfactory describe the observa- 
tions. 

We conclude that all the isochrones have problems with 
masses around 0.6 m©, which is at the transition from p re-MS 
to MS at the age of the Hyades (see, e.g lSiess et"alll200()h 

It is a mass-Mi/ relation for masses above 1 mg, and a mass- 
Mk^ relation for masses lower than 1 m©. 



5. Spatial structure 

The analysis of the spatial structure of the cluster does not only 
require the three-dimensional distribution of the stars as parti- 
cles, but also their individual masses, v 

For the 724 candidates we determined individual masses 
using the following mass-to-luminosity relations: for systems 
with My < 5.38 corresponding to m > mg we take the Mass- 
My relationship for the Hyades derived by [Pinsonneault et al.1 
(2004). For > 5.38 masses are determined via the Mass-M^ 



relationship from the DARTMOUTH isochrones ( Dotter et al 



[2008h . At the transition, we find an insignificant discontinuity 
in the mass determination of less than 0.005 hiq. Finally for 
Mj^. > 5.5, masses are determined via the BCAH isochrones 
(iBarafi^e et aI.LIl99^ . The transition at Mk, - 5.5 creates an in- 
significant step in the estimated masses of 0.01 m^. In the range 
4.0 < < 6.0 masses determined by either of the isochrones 
(DARTMOUTH or BCAH) do not differ by more than 0.03 mo. 
In our application of these mass-to-luminosity relations we ig- 
nore the possible binary nature of the stars in our sample, since 
we have information on binarity only for a minor portion of our 
sample. 

Using the individual masses of stellar systems, we derive 
the cumulative mass function M(rc) depending on the distance 
rc from the cluster centre. This function is shown in Fig. |6] It 
is corrected for contamination as described in Sec. 14.3.11 In the 
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galactic disk the sphere of influence of a gravitational body is 
given by 



XL = 





' GMc "1 


UA(A-B)j 





(3) 



where Xl is the distance of the Lagrangian points from the cen- 
tre, Mc is the total mass inside a distance Xl from the centre, G 
= 4.3 X lO"^pc/m0(km/s)^ is the gravitational constant, A and B 
are Oort's constants, Oq the angular velocity and k the epicyclic 
frequency at the position of the Sun. Here w e use A (14.5 
kmsj'kpc ') and B (-13.0 kms 'kpc ') from iPiskunov et al.l 
(12006) . The distance Xl is often referred to as the tidal radius rt 
of a cluster and we use this definition in the following also. The 
tidal radius rt separates, in general, stars gravitationally bound to 
a cluster from those that are unbound. However, this definition 
is not to be understood that each star inside rt is bound, and all 
stars outside are unbound. Individual stars can, every now and 
then , change from a bound to an unbound state and vice versa 
(see lJust et all l2009l for the fraction of potential escapers with 
Jacobi energy above the critical value in the cluster). It is in- 
herent in the convergent point method, that we find, in the first 
place, stars that are co-moving with the adopted space motion of 
the centre of the Hyades cluster. So, we cannot separate, from 
the beginning, the gravitationally bound Hyades cluster and its 
surroundings. 




Distonce from centre, pc 

Fig. 6. The cumulated mass as a function of the distance from 
the cluster centre rc. The thick black curve is the observed cu- 
mulated mass function as obtained from adding up all the indi- 
vidual masses within a distance rc. The (thinner) red curve is the 
relation from Eq.|3]between the tidal radius and the tidal mass in 
the gravitational field of the Galaxy near the position of the Sun. 



If we would add up the masses of all our Hyades co-moving 
stars within 30 pc from the centre we would get a tidal radius rt 
of 10.5 pc from Eq.[3] Now, it is clear that stars outside 10.5 pc 
are not gravitationally bound, and cannot self-consistently con- 
tribute to the determination of rt itself. In general, the tidal radius 
is the distance rt where M(rt) calculated via Eq.[3]equals the inte- 
grated mass of all stars up to the distance rt from the centre. The 
function shown as the thin red curve in Fig.|6]represents the rela- 
tion from Eq.[3] and the black curve gives the counted integrated 
mass of our stars as a function of the distance from the centre. 
From Fig. |6] we find that the tidal radius comes out to be 9.0 pc 
and the corresponding tidal mass is 276 mg. Uncertainty in the 
determination of rt may come from the fact that we use system 
masses in the case of unresolved binaries, that we neglect the 
contribution from white dwarfs, and that our sample is getting 
incomplete at the lowest masses. The additional mass in white 



dwarfs, found so far in the Hyades , is ver y low. It has been esti- 
mated to be 6.4 m0 by Von Hippell (119981) . We show in Sect. 16.21 
that incompleteness at the low-mass end of the mass function has 
only negligible influence on the tidal radius. We would miss less 
than 8% of the total mass in the unrealistic assumption that the 
logarithmic mass function were constant from 0.25 m© down to 
0.01 m©. The contribution due to unresolved binaries to the clus- 
ter mass is m ore substantial. Based on their radial velocity study, 
iGunn et akl (0988) argue that about 25% of the stars they mea- 
sured are binaries. Assuming the average mass of the secondary 
of 2/3 of the primary's mass. lGunn et al.l (^988) estimate an in- 
crease in the total mass by 17%. Applying this correction would 
give us an increase in the tidal radius of only 5%, from 9.0 pc 
to 9.45 pc. Even a 50% increase in the total mass would only 
enlarge the tidal radius to 10.4 pc. 

Another source of uncertainty for the determination of the 
tidal radius are the adopted values of Oort's constants to char- 
acterise the tidal field of the Galaxy. Using different values for 
the parameters of galactic rotation, we could determine a tidal 
radius between 8.0 pc (using D.o = 30.6, and k = 39.0, both in 
kms ~'kpc~') and 9 . 3 pc ( using A = -B = 13.75 kms^'kpc"') . 
The lPiskunov et al] (l2006h values for A and B give a tidal radius 
in between these two approaches. So, for the rest of the paper, 
we use a tidal radius of rt = 9 pc. In models of star clusters, es- 
pecially in A^-body calculations, the Lagrangian radius for 50% 
of the cluster mass, plays an important role and is called half- 
mass radius i\. Estimated from a tidal mass of 276 m©, the half- 
mass radius is rh =4.1 pc at the present state of evolution of the 
Hyades. 




Fig. 7. The dependence of number density, mass density, and av- 
erage mass per star as a function of the logarithm of the distance 
from the centre rc. From top to bottom: a) The number density 
D of stars in the cluster obtained with a binning of 2 pc in steps 
of Ipc. Outside 10 pc the bin size is increased to get a better 
signal-to-noise ratio. The error bars are Poisson errors, b) The 
observed mass density p (black dots) using the same binning and 
a density distribution (red crosses) from a Plummer model (Eq. 
|4]with M, = 276 m© and rco = 3.10 pc). c) the average mass per 
star obtained as p/D. Its decrease from the centre outbound char- 
acterises the mass segregation in the Hyades, at present time. 
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In Fig. |2]we show the number density D, the mass density p 
, and the average mass per star as a function of the distance from 
the centre rc. The bottom part of Fig.|7]reveals the mass segrega- 
tion in the Hyades; the average mass per star decreases from 1 .4 
mo in the centre to about 0.5 m© at the tidal radius rt. We discuss 
this in more detail in Sect. 16.21 together with the description of 
the mass function. The mass density of the Hyades, shown in the 
middle part of Fig. |7] is fitted to a Plummer model (.Plummen 
[Ml, where the mass density follows the equation 



Table 2. Shape parameters of the observed Hyades O, and the 
model cluster M. See a detailed description in the text. 



P(rc) 



3M, 



1 



47rreo^^ [l + (rc/rco)2]5/2' 



(4) 



where rco is the so-called core radius of a cluster Using the tidal 
mass of M, = 276 m© the best fit to the observed density distri- 
bution (black dots) in Fig.Qis obtained with a core radius of rco 
= 3.10 pc. The core radius is easily seen in the two upper panels 
of Fig.Qas the distance where the slope in the density and mass 
distributions significantly changes. The corresponding Plummer 
model (crosses in Fig.|7]i shows excellent agreement with the ob- 
servations. The model has a central mass density of 2.21 mopc^^* 
which coincides well with our innermost point of 2.14 mopc"^. 
The ratio of the half-mass radius rh to rco in the Hyades is 1.32 
which also is in remarkable coincidenc e with the theoretical ra - 
tio of 1.3048 fo r a Plu mmer model (see lSpitzer & Shulll [19751) . 

iGunn et aH (Il988h already found that the cluster conforms 
well to a Plummer model. Taking into account the incomplete- 
ness of their sample, they argue for a total mass of 390 m©. They 
use the projected density distribution to determine a core radius 
rco = 3.15 pc (0.07 radians) and a central density of 2.97 mopc"^, 
the latter being about 30% higher than ours. As the core radii rco 
in Gunn's model and ours practically coincide, the difference in 
centra l density comes sol ely via the tidal mass of the cluster So, 
either iGunn et al] (1 19881) have overestimated the tidal mass of 
390 mo, or our stellar masses are underestimated essentially due 
to binarity. Using Gunn's estimate on binarity (25%) we could 
account for an increase of 17% in total mass to get 323 m©. A 
mass of 390 m© would require an increase of 37% in mass for 
all our stellar systems (more than 50% bina r ies), w hich, how- 
ever, could also be possible. iPerrvman et al.l (11998) also fit the 
observed density distribution to a Plummer model and find a 
smaller core radius of 2.9 pc and a central density of 1 .8 mopc"^. 
This is not surprising as the Hipparcos sample is incomplete al- 
ready at 0.6 to 0.7 mo. 

In the following we call the region inside rco = 3.1 pc the 
core, the region between rco and rt the corona. In a volume be- 
tween 1 and 2 rt, called halo, we find another 100 mo in ob- 
jects co-moving with the Hyades, but probably not gravitation- 
ally bound. Outside the halo up to the distance of 30 pc from the 
cluster centre there is another 60 mo in co-moving stars. 



5.1. The shape of the cluster 

Since the convergent point method provides an individual dis- 
tance from the Sun for each Hyades star, we can study the three- 
dimensional distribution of cluster members in more detail. The 
distribution of the 724 Hyades candidates in galactic rectangular 
coordinates is shown in Fig. l8^. 

For comparison, we show in Fig. ISJ) a model distribution 
of st ars predicted by the A^-body simulation (iKharchenko et al.L 
12009 ) of an open cluster originally spheroidal, with an initial 
mass of 1000 mo and the Salpeter IMF down to 0.1 m©. The 
model cluster is moving on a circular orbit in the external tidal 
field of the Milky Way. It is 650 Myr old and located close to 
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[0.70; 0.72] 


[0.43; 0.45] 


[-i;2] 


[25; 36] 


[2; 7] 


0/18 
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the observed location of the Hyades (fo r more details of the 
model parameters, see IKharchenko et al.L 12009 ). It is appropri- 
ate to note that the model has not been specifically tailored to fit 
the Hyades, but it is a general model for the evolution of a pro- 
totype open cluster in the disk at 8.5 kpc. The model does not 
take into account that the Hyades could undergo some irregular 
events (e.g., encounters with molecular clouds) during their life. 
As far as the observations are concerned, we note that our sam- 
ple can be contaminated beyond rc = 9 pc. However, the model 
and observations show quite a good agreement in general. The 
main feature, an elongated shape with the minor axis directed 
towards the Galactic Poles, is clearly seen in both distributions. 

In the following we analyse and compare the shapes of the 
observed Hyades cluster and the model cluster The principal 
axes of the cluster ellipsoid are given as the eigenvectors and 
eigenvalues of the matrix formed by the 2nd order momenta 
Mf^ , ■ • ■ M., of the space coordinates of the cluster stars 



M 



M,, M,,. M,, \ 

M,,, M„, My, 

Mf. My, M,^ 



The coordinates x, y, z are counted relative to the cluster centre in 
the galactic coordinate system as defined in Sect.|3] Let A, B, C 
be the eigenvalues of M with A being the highest, C the lowest 
eigenvalue. The corresponding eigenvectors are denoted by 

b and"c^. 

The results are given in Table 13] In the first column we in- 
dicate the data set used for the calculations (Observations or 
Model, within 30 pc, 18 pc, or 9 pc). Columns 2 and 3 give 
the axis ratios. Columns 4, 5, and 6 describe the orientation of 
the ellipsoids in the XFZ-system where if/ is the angle between 
the X-axis and the projection of the semi-major axis ~i on the 
XZ-plane, ip is the angle between the X-axis and the projection 
of ~d on the XF-plane, and 9 is the angle between the F-axis 
and the projection of b on the KZ-plane. The numbers in brack- 
ets below each parameter give a measure of confidence. These 
are the lowest and highest values obtained for a given parame- 
ter from four subsets of data in the corresponding range of rc. 
The first two subsets were extracted from the sample arranged 
with increasing distance of stars from the cluster centre by sort- 
ing into two groups of even and odd entries. These two subsets 
contain, therefore, half of stars of the original sample and are 
supplementing each other The other two subsets were extracted 
in the same way but from the sample arranged with increasing 
distance of stars from the Sun. 
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Fig. 8. Observed (a) and modelled (b), (c) spatial distribution in the Hyades. Blue points are for stars with masses m > 0.5mo, red 
points for m < O-Sm©. In the modelled distribution, the dynamically unbound particles are marked by black crosses (b) or black 
points (c). Two large black circles in (a) and (b) indicate the distances r^ from the cluster centre of 9 pc and 18 pc, whereas in (c) 
the circles correspond to rc = 9 pc, r^ = 18 pc, and rc - 30 pc. The last one restricts the area considered in this study. 



Although the model particles are distributed more regularly 
than the real stars, both distributions show similar parameters 
and tendencies. If the area within rc - 30 pc is considered, the 
semi-major axis A of the observed as well as of the modelled 
ellipsoids is about a factor of two larger than the semi-minor 
axis C (C/A = 0.58 for O/30, and C/A = 0.44 for M/30). The 
axis is located in the XF-plane (if/ a 0°) but inclined to the 

X-axis by an angle (p of about 30°, whereas the second axis b 
seems to form a small angle with respect to the XF-plane. With 
decreasing rc the distribution of stars becomes more spherical, 
and the orientation angles lose their meaning. 

An interesting feature in Fig. |8];, the tails parallel to 
the F-axis, is predicted by the model (see also Fig. 3 in 
iKharchenko et al.L l2009l) . These tails spread out to ^ 700 pc 
from the cluster centre and they are sparsely populated by for- 
mer cluster members, mainly low-mass stars with space veloc- 
ities pretty different to that of the cluster itself. Therefore, even 
if the tails would not be destroyed by events not taken into ac- 
count by the model, they are difficult to segregate from the field 
by use of kinematic criteria alone. From this point of view, all- 
sky surveys of chemical compositions of stars would be more 
promising to detect former members of an individual cluster. 



6. The mass and luminosity functions 

In this paper we derive the masses of our Hyades candidates 
(stellar systems) from a mass-luminosity relation (MLR). Hence 
the luminosity function and the mass function are not indepen- 
dent of each other, and cannot be discussed independently. From 
the point of view of an observer, the luminosities are the pri- 
mary observable quantities and the masses are derived thereof. 
However, from theory, the primary astrophysical parameter is the 
initial mass function (IMF) which evolves into the present-day 
mass function (PDMF). In consequence, we shall epitomise the 
observed luminosity function in the next subsection, and discuss 
the PDMF in more detail. 



In section |5] we described the mass-luminosity relation we 
adopted to convert absolute magnitudes into masses. It is a mass- 
Mv relation for masses above 1 m^, and a mass-MK, relation 
for masses lower than 1 m©. Our derived luminosity and mass 
functions in this section were corrected for contamination as de- 
scribed in Sect. gXU (7.5% for 9 pc < rc < 18 pc and 30% at 
18 pc < rc < 30 pc). In section l4r2l we found that our sample is 
getting incomplete between Mk, =7.3 and Mk^ = 8. Translated 
into mass, incompleteness in the mass functions sets on some- 
where between 0.25 and 0.17 m©. The onset of incompleteness 
is shown by the dashed vertical lines in Figs. l9l and [TOl 

6. 1 . The luminosity function 

The K-band luminosity function (KLF) is defined as the number 
of stars per magnitude bin and volume element 



dN(MK,) 
dMK ' 



(5) 



The KLF is mostly used to characterise the stellar population in 
young or even embedded clusters. On the other hand, JahreiB 
(priv. comm.) has recently determined the KLF of the stars in 
the solar neighbourhood from the CNS4. We compare the KLF 
from the Hyades both with that of young clusters and with the 
older local population. 

In different regimes of magnitude, the KLF is usually char- 
acterised by a power law in the form 



dMv 



(6) 



with a a s the slope of the power law. For you ng open clusters 
(see, e.g. iDevine et al.L l2008t iJose et all l2008l as an overview) 
the slopes a of the KLF lie between 0.3 and 0.4 in ranges of Mk, 
between and 3, whereas the slope of the KLF of CNS4 (see Fig. 
|3]l is steeper, about 0.5, in the corresponding magnitude range. 
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To be independent of binning we construct the luminosity 
function in steps of 0.1 mag with a binsize (kernel) of 1.5 mag. 
We also distinguish different regions of the cluster as defined in 
Sec. |5] the core region up to rc = 3.1 pc from the centre, the 
corona between core and tidal radius of 9 pc, the halo in the 1 
to 2 tidal radii shell, and the co-movers outside 18 pc up to 30 
pc from the cluster centre. We show the KLF for the different 
regions of the Hyades in Fig.|9] 



D 

E 



-2 - 



I 



u 

CL 



-3 



-4 




Fig. 9. The luminosity function of the Hyades in different re- 
gions. This luminosity function is constructed as a moving aver- 
age with a rectangular kernel of 1 .5 magnitudes in steps of 0. 1 
mag. The upper (blue) curve refers to the core region inside 3.1 
pc from the centre, the second (cyan) curve to the corona (3.1 
pc to 9 pc), the third (green) to the halo (9 pc to 18 pc), and the 
fourth (red) to the co-movers outside 18 pc. The dashed verti- 
cal line at Mk, - 7.3 marks the onset of incompleteness of our 
survey. The units are stars (systems) perpc"-' and 1 mag. 



Looking at the KLF of the different regions one feature im- 
mediately catches the eye: the luminosity function in the core of 
the Hyades is strikingly different from the other regions and also 
from the field star LF in Fig. |3] We observe a strong excess of 
luminous stars with Mk, < 5. This is a manifestation of mass 
segregation already seen in Fig. Q 

The maxima of the PDLFs in the different regions are 
reached at about Mfc^ - 7, except for the additional maxima in 
the core at Mk^ = 4 and Mfc^ = 1.5, caused by mass segregation. 
The local KLF from Fig. [3] is still increasing at Mk^ = 7, and 
reaches its maximum at Mk^ = 9 or later This implies a lack 
of low-luminous stars in the Hyades compared to the field with 
all caution, because incompleteness in our sample starts at My . 
= 7.3. In the optical regime, according to iKroupa et al.l (Il990h . 
the maximum of the field star LF is located at My =12 which 
corresponds to Mfc^ - 7.3. This is at a brighter magnitude than 
the JahreiB (priv. comm.) KLF. We have no explanation for this 
difference in the field star LF. 



At absolute magnitudes 3 < M[c, < 5 we find local maxima 
in all KLFs shown in Fig.|9l In the KLF of CNS4 we do not see 
a maximum but a flattening of the slope to o- = 0. 1 . In the opti- 
cal th i s phenomenon is kno wn as the Wielen dip dKroupa et all 
Il990h . IKroupa et al.1 (Il990h located it at m = 0.6 mo. They ex- 
plained the dip as a consequence of an inflection of the mass- 
luminosity relation which appears due to the formation of nega- 
tive H-ions impacting the atmospheric opacity. A mass of 0.6 m© 
corresponds to Mk^ = 5.1 in the MLR from Sect. |5] We observe 
strong Wielen dips (even local minima) in the luminosity func- 
tions in all regions of the Hyades in Fig.|9] strongest in the core 
region, where mass segregation is responsible for an apparently 
stronger Wielen dip. 

For 1 < Mk, < 4, we find the slope of the luminosity function 
in the halo to be a = 0.26, i.e marginally flatter than in young 
open clusters and much flatter than in the field. This behaviour in 
a zone of co-moving, gravitationally unbound stars, shows that 
we find there an ensemble of "former" members representing 
the IMF of open clusters at least in this magnitude range. On 
the other hand, the slope is remarkably different from that in 
the field star KLF meaning that the halo is exclusively related 
to the cluster and not to the field. This can be considered as an 
alternative confirmation of low field star contamination (cf . Sect. 

6.2. The mass function 

The mass function f (m) of a cluster is defined as the number of 
stars within a mass interval between m and m + dm 



^(m) 



dN(m) 
dm 



(7) 



Alternatively, we can describe it as a function of the logarithm 
of mass 



dN(lo2 m) 

^(log m) = \ ^ ' = m X ^(m) X In(lO). 
d(logm) 



(8) 



So, ^(log m) gives the number of stars (systems) in the in- 
terval between logm and logm + dlogm. In Figs. [TOl and [TTI 
we chose d log m to be 0. 1 . The different regimes in the mass 
function ^(m) are usually approximated by a power law with 
an exponent -a which transforms into -a + 1 for the logarith- 
mic mass function ^(log m). So, even if we show the logarithmic 
mass function in the figures below, we always refer to the expo- 
nent -a of ^(m). 

Fig. [To] shows the mass function in different regions of the 
cluster The lowest function (black) refers to the core region 
(rc <3.1 pc). The mass function peaks sharply at about 1 m©. 
Towards larger masses the MF drops as a power law with expo- 
nent a = 2.7, even steeper than Salpeter The most striking effect 
is the steep decline {a = -0.05) towards lower mass stars, again 
a clear indication for mass segregation. 

The mass function for all stars within the tidal radius of the 
cluster re < 9 pc has its less prominent maximum at 0.9 m^, and 
follows about the same power law {a - 2.7) towards higher mass 
systems as in the core (provided that we neglect the four giants in 
the most massive bin). Towards low-mass stars the logarithmic 
mass function is almost flat (a = 0.93) down to the complete- 
ness limit. This mass function refers to the gravitationally bound 
stars of the Hyades cluster Adding the formerly bound members 
with rc > 9 pc shifts the maximum further towards lower mass 
stars. However, the logarithmic mass function for all 724 stars 
is almost flat in the mass range from 1 m© to 0.2 m©. The mass 
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Fig. 10. The present-day mass function (PDMF) of the Hyades. 
The lowest histogram (filled black dots) shows the PDMF in the 
core of the cluster (rc <3.1 pc), the magenta triangles mark the 
PDMF within the tidal radius (rc <9 pc), the blue squares stand 
for the PDMF in the region rc < 1 8 pc, and finally the filled red 
squares are the PDMF of all 724 stars within 30 pc from the cen- 
tre. The PDMFs are corrected for field star contamination. The 
histograms give the actual star counts per logarithmic mass bin 
of 0.1. The dashed vertical line at about 0.25 niQ marks the on- 
s et of incomplet eness of our survey. For comparison, we show 
a iKroupal (l200ll) mass function (solid line) with a - 2.35 for 
masses larger than 0.5 m©, and a = 1.35 below 0.5 m©. As an 
other extreme we show a mass function with slope a - 3.05 
(dashed line) fitted to the high-mass portion of the observed 
PDMF. 



should point out again that in the mass functions presented here 
binaries are not resolved, and the masses are based on a specific 
mass-luminosity relation used to convert absolute magnitudes to 
masses. This makes us cautious, not to interpret the slope of the 
mass function at high masses in more depth. 



1000 
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Fig. 11. The mass function of the Hyades from this paper (black 
histogram with filled black dots). The blue curve (open squares) 
shows the Hyades PDMF from Bouvieretal. (2008), the red 
histogram wit h triangles is the PD MF from the 120 Myr old 
Pleiades from lMoraux et ak ] (I2003 h. and finally the green curve 
(filled squares) is the PDMF of Praesepe from Boudreau lt et al.l 
(2010). The thin solid line shows again a Kroupai (1200 iF mass 
function. 



function which we determined is incomplete for masses lower 
than about 0.25 m^. This has, however, not much influence on 
the total mass of the cluster If we assume that the logarithmic 
mass function were constant below 0.25 m© down to 0.01 m©, 
this would add merely 8% to the total mass of the cluster 

In FigfTOl we compare the observed mass functions with 
a IKroupal T2OO lb IMF shown as the solid line. This IMF has 
a = 2.35 for masses larger than 0.5 m^, and a = 1.35 below 
0.5 mo, down to the brown dwarf limit of 0.08 m©. Using the 
Kroupa IMF as a tangent to the observed mass function, we can 
estimate a minimum initial mass of the cluster, provided that all 
originally formed stars of about solar-mass are still contained 
within a radius of 30 pc. With these assumptions we get an 
initial cluster mass of 1100 hiq in 2300 stars (stellar systems). 
From this simple calculation we can conclude that about 40% 
of the original mass of the Hyades is still contained in the im- 
mediate neighbourhood of 30 pc around the centre. The remain- 
ing mass has left f urther into tails of the form predicted by, e.g. 
iKharchenko et al.l (i2009) . More refined statements on the initial 
parameters of the Hyades are expected from dedicated A^-body 
simulations which are underway. 

A better fit to the high-mass end of the observed mass func- 
tion requires a steeper slope of a = 3.05 in a mass range from 
1 to 2.5 m©. This is shown as the dashed line in Fig. [TOl We 



In Fig.[TT|we compare our PDMF of the Hyades with earlier 
results. This mass function (black histogram with filled black 
dots) is obtained within the 30 pc volu me around the centre 
(top curve in Fig. \M- The PDMF from Bouvier et all (l2008l) 
is shown as the blue histogram (open squares). On the high mass 
part it is derived fro m the Prosser & Stauffer (P&S) database. 
iBouvier et al ] (120081) argue that the P&S database is essentially 
complete down to 0.1 m©, but having a 27% contamination by 
field dwarfs below 0.3 m©. The slope of their mass function be- 
tween 0.2 and 1 m© is found to be a = 0.6. Compared with our 
results this indicates that incompleteness in the P&S database 
begins near 1 m© and increases towards lower masses. We at- 
tribute this to the fact that the P&S database is not based on a 
deep all-sky survey. As fainter stars show less concentration to 
the centre (mass s egregation) the y may be missed if the field-of- 
view is restricted. iBouvier et al.l (l2008) detected 3 new very low- 
mass stellar members, and constructed a PDMF down to 0.08 
m0. Nevertheless, their PDMF practically shows a cut-off' at 0.1 
m0. We neither can confirm nor reject this on the basis of our 
analysis using the CU catalogue because of incompleteness be- 
low about 0.25 mg. Work is going on to exploit the faint PPMXL 
stars not in the CU subset to reach an all-sky completeness of the 
Hyades stars at least down to 0.1 niQ. Hogan et al. (2008) did not 
try to derive a mass function from the 12 L-dwarfs in their sur- 
veyed 275 deg- field. If we apply our Mass-M/j-, relationship, the 
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12 candidates range from 0.078 down to 0.06 hiq. A larger field 
is needed to decide if the presently observed strong cut-off at 0. 1 
m© will be confirmed or not. 

Fig. [TT] also shows the PDMF of the much younger 120 
Myr) Pleiades cluster from iMoraux et all (120031) . The maximum 
of the Pleiades PDMF in the logarithmic mass function is be- 
tween 0.3 and 0.4 m© compared to about 0.7 to 0.8 m© in our 
mass function of the Hyades. If both clusters had a similar IMF 
then mass loss works more effectively the less massive the stars 
are. The similarity of the mass functions in the region above 1 
mo suggests that this effect is less prominent for the more mas- 
sive stars in a cluster A further similarity in the Hyades and 
Pleiades mass functions is the decline below 0.1 m©. 

The green curve (filled squares) is the mass function of 
Praesepe recently derived bv lBoudreault et aD(l2010l) . Given that 
Praesepe is practically of the same age as the Hyades, this mass 
function is quite unusual. With an a = 1.9 between 0.2 and 
0.7 m0, the contribution of low-mass stars to the total mass of 
Praesepe is much more important than in the case of the Hyades 
and even the Pleiades. Contrary to the latter, it would mean that 
Praesepe could hold its low-mass stars in bound stage. 

T he thin s olid l ine in Fig. [TT] is again the universal IMF 
from iKroupal (1200 lb . Prima facie one gets the impression that 
an IMF like Kroupa's evolves via a Pleiades PDMF (120 Myr) 
to a Hyades PDMF due to dynamical and stellar evolution. 

7. The internal velocity dispersion 

For a cluster in a statistically steady state (IChandrasekharl 0*9421) 
the virial theorem 2T = -f2 holds (with T being the total kinetic 
and Q. the total potential energy of the system). For compact, iso- 
lated, massive and long living systems like globular clusters the 
assumption of statistically steady state is, in general, justified. 
However, this is not trivial if we were to consider a relatively 
old, low-mass open cluster dipped into the Galactic disk. 

Strictly speaking, the virial theorem only holds for the 
cluster as an ensemble. For each well-defined subset of stars 

{nil, ,ms] we can postulate that the velocity dispersion is 

given by 

<v^>=0.5g2^/Z-' (9) 

where m, is the mass of an individual star in the subset, and nij 
runs along all stars within the tidal radius of the cluster, r,j is 
the distance between stars i and j, G is the gravitational constant, 
and < v^ > is defined as 

< v^ > x^mi = ^miVi^ (10) 

The summation in Eq.|9]extends over all pairs of stars with i j. 
If the subset consists of all stars within the tidal radius then Eq. 
Ogives 2T = -Q. 

We use here the actual three-dimensional distribution of the 
Hyades stars and do not approximate the cluster by a Plummer 
model, albeit the fit of the Hyades to a Plummer model is excel- 
lent (see Sect|5]l. The assumption of a statistically steady state, 
if at all, is only justified within the tidal radius of the cluster. We 
applied Eq. |9]to two distinct subsets of the Hyades, one contain- 
ing all stars within the core radius rco (rc < 3.1 pc) and the other 
within the tidal radius rt (rc < 9.0 pc). The three-dimensional 
velocity dispersion comes out as 0.42 kms"' and 0.32 kms"', re- 
spectively. Putting an artificial particle in the centre we obtain 
the central velocity dispersion < v^ >'^^ = 0.45 kms"'. 



iGunn et aP (Il988l) used their fit of the Hyades cluster to a 
Plummer model (cf. Sect|5]l to get 0.5 1 kms ' for the central, and 
0.40 kms' for the average velocity dispersion (or0.23kms"' for 
1-D). The higher values can be explain ed by their assumption 
of a higher (tidal) mass for the Hyades. iPerrvman et al.l (Il998h 
derive a central velocity dispersion of 0.36 kms"' from their fit 
to a Plummer model, which is too low because of their lower 
mass estimate for the cluster 

In principle, an empirical determination of the velocity dis- 
persion is possible by exploiting the observed velocity distribu- 
tion of cluster members. However, during their life, open clus- 
ters are permanently losing members part of which remain in the 
vicinity of their parent cluster and form sparse haloes and tails 
of co-moving stars. This makes it difficult to separate "actual" 
members by use of velocity observations and properly determine 
< v^ >'''^. For open clusters, one expects velocity dispersions 
ranging from a few hundred meters per second to 1-2 kms ' i.e., 
at a level of accuracy of radial velocity measurements and/or tan- 
gential velocities of nearby stars. Empirically determined < v^ > 
are thus sensitive to the adopted mean errors of input data used 
for computing int ernal spatial velocities of cluster members (see 
e.g.. Appendix in lPiskunov et ali|2007h . 

In the past, several attempts were undertaken to derive the 
velocity dispersi on of th e Hyad es from proper motions and 
radial velocities. ISchwanI d 19901) determined the velocity dis- 
persion from 44 FK5 stars. Using Vj^, he obtained an Id- 
velocity dispersion of 0.71 + 0.41 kms Based on only the v^- 
component, the velocity dispersion comes out to be 0.3 kms ' 
which corr esponds t o about 0.5 kms ' for the three-dimensional 
dispersion. ISchwanI ((1990) explained this considerable differ- 
ence by larger uncertai nties in v^-component of the velocity data. 
iPerrvman et al.l (Il998h used the Hipparcos parallaxes and proper 
motions, but were unable to determine the velocity dispersion, 
because the accuracy of the measured trigonometric parallaxes 
resulted in a too large mean error of the tange ntial velocities. 
'de Bruiine et al.l(l200lh used the same data set as IPerrvman et al.l 
(1998^ but took the formally better secular parallaxes instead of 
the trigonometric parallaxes. They achieved a Id-velocity dis- 
persion of 0.3 kms"' after consecutive rejection of stars with 
large offsets in the observed vs. computed tangential velocities. 
In the following we estimate the three-dimensional velocity dis- 
persion in the Hyades via the one-dimensional velocity disper- 
sion given by the velocity component perpendicular to the direc- 
tion to the convergent point v^. The distribution of v^ versus the 
distance of stars from the cluster centre is shown in Fig. [12] Stars 
in red circles are outliers and are discussed a few paragraphs be- 
low. 
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Fig. 12. The distribution of the tangential velocity components 
Vj. perpendicular to the direction to the convergent point as a 
function of the distance from the cluster centre. The nine stars 
in red circles are outliers which are not considered in the deter- 
mination of the velocity dispersion in Table [3] The vertical Unes 
mark the core (dashed) and tidal (solid) radii. 
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Let V, be the residual space velocity of the i-th member with 
respect to the bulk motion of the cluster, and cr^ its individ- 
ual mean error We define the square of the velocity dispersion 



< V > as 



2 

< V > = 



(11) 



where n is the number of stars in a given region. Provided that 
the space velocities of the cluster stars are randomly oriented, 
the first term on the right-hand side can be computed from the 
velocity components Vj. and we can re-write Eq.fTTIas 



vj =< v^ >= 3 X (< vi > - < a-l >). 



(12) 



The mean errors ctvj (iTrumpler & WeaveR Il953l) are computed 



Table 3. Three-dimensional velocity dispersions derived from 
the potential energy (Eq. |9]l as well as computed from the ob- 
served velocities of stars in 3 mass ranges: I (m > 1.45 m©), II 
(1.05mo<m< 1.45 mo). III (0.71mo<m< 1.05 mo). 





theoretical 


I 


II 


III 




<v>'/2 


< V>'/2 


<V>'/2 


<V>'/2 


pc 


kms"' 


kms"' 


kms"' 


kms"' 





0.45 








< 3.1 


0.42 


0.81 ±0.13 


0.76 ±0.15 


0.74 ± 0.21 


< 9.0 


0.36 


0.77 ±0.10 


0.88 ±0.11 


0.82 ±0.13 


> 9.0 






2.68 ± 0.38 


2.38 ± 0.24 



Vrf(2n)i/2- 



(13) 



Here we are obliged to rely on the quoted cr, and to assume their 
uncertainties to be negligible. 

Since the kinematic candidates of the Hyades have different 
observing history, their individual mean errors crj_ vary signifi- 
cantly, from ^0.1 kms ' for bright (and high mass) Hipparcos 
stars to «; 1 kms"' for the faintest (and lowest mass) stars in our 
sample. Therefore, they must be considered in different bins con- 
taining a sufficient number of stars of homogeneous accuracy. As 
a compromise, we divided our sample into four groups, accord- 
ing to the masses of stars included: masses higher than 1.45 mo, 
between 1.45 m© and 1.05 m©, between 1.05 m© and 0.71 mo, 
and less than 0.71 mo. We had to exclude the last group from 
consideration: there is no chance to use them for the determina- 
tion of the velocity dispersion due to the low accuracy of their 
kinematic data. Further we rejected 5 stars from the calculation 
of < cr^ > because their reported mean errors were significantly 
higher than the corresponding group average. 

For each mass group, the mean values < v^ > were com- 
puted in two ranges of the distance from the cluster centre: 
rc < 3.1 pc, and rc < 9 pc. It is self-evident that stars with high 
I v^ I can cause an oversize increasing of the velocity disper- 
sion. Therefore, in each bin we applied a statistical test to reveal 
outliers at the 95% confidence level. In total, 9 stars were ex- 
cluded from the calculation of the velocity dispersions. They are 
marked by red circles in Fig. [12] All these stars are doub l e stars . 
Of these, six are found as A/i-binaries in Wi elen et al.l (U999I). 
The binary nature of the AyU-stars is revealed from the signifi- 
cant diff'erence between their short-scale (Hipparcos) and their 
long-scale (FK5, Tvcho-2) proper motions. One more Aju-binary 
is found in iFrankowski et al] (|2007|) . Another star is marked as 
double in the Hipparcos Double and Multiple Systems Annex, 
and the last one is WDS 04285-1-174 2 from the 201 1 ver sion of 
the Washington Double Star catalog ( iMason et "all 12001 1) . 

For comparison, we also computed the velocity dispersion 
of stars at larger distances from the cluster centre by applying 
the same approach. Since there are only 8 stars at rc >9 pc with 
masses larger than 1.45 m©, the calculations were carried out 
only for two mass groups. The results are given in Table |3] to- 
gether with the theoretical velocity dispersion computed above 
from the virial theorem. 

The increase of the velocity dispersion outside 9 - 11 pc 
from the cluster centre can be akeady observed in Fig. [12] and 
Table [3] confirms this impression. Partly, this may be explained 
by an increasing contamination by field stars. But more proba- 
bly, this finding indicates the limiting distance for the gravita- 
tionally bound part of the Hyades cluster The distance 9 - 1 1 pc 



is in good agreement with the conclusion on the cluster size we 
obtained from the density and mass distribution of stars in Sec. 

m 

Within the first two bins (r^ < 9 pc), the velocity disper- 
sion achieves 0.8 kms"' which is higher by almost a factor 
of two than one would expect from Eq. [9] So, a factor of 4 in 
mass would be required to balance this discrepancy. Assuming 
additional mass hidden in double stars, white dwarfs, faint stars 
below the completeness limit of our survey, the current estimate 
of cluster mass of approximately 300 m© can be hardly increased 
actually by more than 50%, a factor of 4 is excluded. 

On the other hand, we can assume that < v^ > is overes- 
timated and/or < o-]_ > is underestimated in Eq. [12] A formal 
calculation shows that a decrease of the observed velocity dis- 
persion by a factor of two can be achieved for the most massive 
stars (group I) by increasing the rms errors of proper motions 
by a factor of 3. However, this seems to be far from real- 
ity, especially for this mass group containing mainly stars with 
Hipparcos data. 

Alternatively, an additional dispersion can be introduced by 
the presence of double stars which can impact the proper motion 
measurements. The effect from neglecting the binary nature of 
the stars is twofold: the velocity dispersion contains a part that 
comes from orbital motion, and the total gravitational energy is 
increased to the higher mass in binaries. Indeed, we found that 
the high velocity dispersion is caused by only 20% of the stars 
with somewhat higher v^ velocities. From a statistical point-of- 
view, however, they cannot be considered as outliers and simply 
be excluded from the calculation of the velocity dispersion. But, 
we note that th i s num ber of 20% coincides well with the results 
by iGunn et ail (Il988l) . who estimated the percentage of double 
stars to be 25% from their measurements of radial velocities of 
stars in the Hyades region and extrapolated this number to the 
complete cluster. As a further possibility to explain the high em- 
pirical velocity dispersion we can assume that a certain number 
of stars is just escaping the cluster Simulations have been started 
to quantitatively investigate the importance of this eff'ect. 

In the discussion above we used the observed distribution of 
stars (with their masses) to derive the allowed velocity disper- 
sion for the cluster to be in equilibrium. That is, the mass of the 
cluster was obtained by counting the members with their individ- 
ual masses. At first sight, see Table [3]the observed velocity dis- 
persion seemed to considerably exceed the allowed one. But, as 
we argued above, there are good reasons that both the observed 
distribution of masses and the observed velocity dispersion are 
consistent with each other. 

Let us, at the end of this chapter, add a remark on the 
estimation of masses from an observed velocity distribution. 
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IChandrasekharl (Il942h introduces the "average" radius R of the 
cluster such that the equation 

M<v2>=-— (14) 

holds. Here M is the total (gravitationally bound) mass of the 
cluster This is fulfilled for 

J_^J_^m^ (15) 

where the summation is now carried out over all m,, nij with / 
j. The so-called "average" radius of the cluster is determined via 
the weighted mean over the inverse distances between individual 
stars. For the Hyades with their tidal mass of 276 we find R 
= 4.53 pc, which in our case happens to be 1.5 times rc„. The 
measured velocity dispersion can be used to determine the mass 
(virial mass) of an open cluster via 

M = I; < v^ > R. (16) 
G 

The two terms on the right-hand side of Eq(T6] < v^ > and R 
depend on the sample for which the velocity dispersion is mea- 
sured. This should be taken into account when deriving the mass 
of a cluster from its velocity dispersion. 

8. Summary 

We used a subset of PPMXL, derived from a combination of 
PPMXL proper motions with UCAC3 and CMC 14 observations, 
and augmented by CMC 14 photometry to perform the most 
complete three dimensional census of the Hyades to date. This 
census yielded 724 stellar systems as candidates moving with the 
same space motion as the cluster These candidates have been 
found from an application of the convergent point method, and 
we verified the resulting secular parallaxes by photometric paral- 
laxes. Although we tested the full catalogue for stars co-moving 
with the Hyades centre, we restricted our analysis to stars with 
derived distances smaller than 30 pc from the cluster centre. We 
did not search for new white dwarfs, as they may be at the de- 
tection limit of the 2MASS, and need a separate approach. Our 
census is complete down to Mk^. ~ 7.3, corresponding to masses 
of 0.25 m0. Incompleteness for more massive stars is no longer a 
question of limited field-size; the only reason for incompleteness 
may be possibly erroneous proper motions and/or photometry. 

The secular parallaxes allow us to construct an empirical 
Colour-Magnitude Diagram (here Mk^ vs. J - Ks) from the 
brightest stars down to Mfc^ ~ 9. Comparing this with 3 differ- 
ent theoretical isochrones relevant to the metallici ty and age of 
the Hy ades we find that only the isochrones from lBaraffe et al.l 
( Il998l) more or less describe the loci of the Hyades stars fainter 
than about Mk - 5 in the CMP, whereas the isochrones from 
iMari^o et al.1 (120081) and jPotter et all (l2008h predict J - K, 
colours which are too blue. 

Using the mass-luminosity relation described in Sect. |5]we 
derive M(r), the cumulative mass as a function of the distance 
from the cluster centre. M(r) is increasing up to our arbitrarily 
set limiting radius of 30 pc. If we consider the tidal field of the 
Galaxy via Eq. [3] we derive a tidal radius of about 9 pc, con- 
sistent with earlier determinations. The M/L relation gives us 
system masses for unresolved binaries. Using these individual 
masses we find that 276 (364 stars) are bound within the 
tidal radius of the cluster, another 100 m© (190 stars) are found 
in the halo up to 2 tidal radii, and another 60 m© (170 stars) are 



found as co-movers outside 2 tidal radii. Here the masses are cor- 
rected for contamination, whereas the number of stars is given as 
counted. 

The spatial distribution of the Hyades has a clear ellipsoidal 
shape considering both bound as well as unbound stars. The 
largest axis is almost parallel to the galactic plane, and forms 
an angle of about 30 deg with the galactocentric radius vector 
(in the direction of the galactic rotation). The ellipticity and the 
orientation of the principal axis may already reflect the onset of 
tidal tails as predicted by the models. 

Within its tidal radius (rt - 9pc) the spatial distribution fol- 
lows excellently a Plummer model with core radius rco = 3.1pc 
and a central density of 2.21 m^pc^^. The half-mass radius of the 
cluster is rj, - 4.1pc, and the ratio rh/rco is 1.32, also in excellent 
agreement with the model value. The "average" cluster radius R 
determined from the distances between the stars is found to be 
4.53 pc. 

Mass segregation is clearly seen in the Hyades. The average 
mass per star drops from 1.1 m© at the centre to 0.5 m© at the 
tidal radius. This mass segregation is also seen in the Ks-hand 
luminosity function KLF. In the core the KLF shows a peak at 
Mk^ ~ 4 which is more prominent than the other at Mfc^ - 7, 
whereas this is vice versa in the corona. We see clear indications 
for the presence of the so-called Wielen dip at about Mk^ - 5. 
Mass segregation exaggerates the Wielen dip in the core, while 
it is present in all other regimes (corona, halo and co-movers) as 
a flattening of the slope of the luminosity function. 

Compared with the KLF of the field stars in the solar neigh- 
bourhood (Fig. O the Hyades KLF is generally flatter than the 
field KLF in all regions of the cluster in the range 1 < Mk, < 
6. The increase in the field KLF for Mk, > 7 is not seen in the 
Hyades KLF with the caution that our sample is becoming in- 
complete there. 

The mass function of all co-moving stars (bound or un- 
bound) is complete down to 0.25 m©. About 94% of the total 
mass is in stars more massive than 0.25 m©. Hypothetically, even 
a constant logarithmic mass function at lower masses between 
0.25 m0 and 0.01 m© (for which we so far see no indication) 
would add only 8% to the total mass we determined. From liter- 
ature there is no indication for a considerable amount of mass in 
white dwarfs. Hence, hidden mass could only be in unresolved 
binaries, and, indeed, we got hints to this from the observed high 
velocity dispersio n. Adopting the estimate of 25% binaries from 
iGunn et an(ll988h . the total mass would increase by 17%. 

Mass segregation is also seen impressively in the mass func- 
tion (Fig. [Tot. Whereas in the core, the PDMF sharply peaks at 
about 1 m0 and steeply drops to both sides, the number of lower 
mass stars increases up to the tidal radius to give a pr actically 
flat l ogarithmic mass function. If the IMF had been a iKroupal 
(12001.) IMF we see strong dynamical evolution on the low-mass 
side. This would hold even more for a Salpeter IMF. A compar- 
ison with the mass function of a much younger cluster like the 
Pleiades supports this conclusion. Adopting a Kroupa IMF we 
derive a minimum mass of 1 100 for the cluster at formation. 

The distribution of the velocity dispersion is clearly different 
inside and outside the tidal radius of the cluster The observed 
velocity dispersion within the tidal radius is not sustained by the 
total gravitational energy of the cluster A relatively high frac- 
tion of binaries is needed to overcome this discrepancy, because 
binaries not only increase the total mass, but their orbital mo- 
tions also may strongly corrupt the proper motions used for the 
determination of the velocity dispersion. Our investigation on 
possible contamination by local field stars has shown that the 
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population outside the tidal radius is strongly related to the clus- 
ter itself and consists of stars that left the cluster in the past. 

Putting all together, we observe our neighbouring open clus- 
ter as an evolving object, which is permanently losing bound 
members that we can pin-point up to a certain distance once they 
left the cluster. The results that we obtained here for the spatial, 
the velocity and the mass distribution of the Hyades give start- 
ing conditions for ongoing A^-body simulations to reveal the past 
and the future of this cluster 
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